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Abstract We performed an extensive non-LTE analysis of the neutral sodium lines 
of Na I 5683/5688, 5890/5896, 6154/6161, and 8183/8195 for disk/halo stars of F- 
K type covering a wide metallicity range (—4 <^ [Fe/H] <^ 0.4), based on our own 
data as well as those collected from the literature. For comparatively metal-rich 
disk stars (—1 <^ [Fe/H] < 0) where the weaker 6154/6161 lines are best abundance 
indicators, we confirmed [Na/Fe] ^ with an "upturn" (i.e., a shallow/broad dip 
around —0.5 <^ [Fc/H] <; 0) as already reported by previous studies. Regarding 
the metal-deficient halo stars, where the much stronger 5890/5896 or 8183/8195 lines 
suffering considerable (negative) non-LTE corrections amounting to 0.5 dex have to be 
used, our analysis suggests mildly "subsolar" [Na/Fe] values down to ^ —0.4 (with a 
somewhat large scatter of ^ ±0.2 dex) on the average at the typical halo metallicity 
of [Fe/H] ^ —2, while they appear to rise again toward a very metal-poor regime 
recovering a near-solar ratio of [Na/Fe] ~ ([Fe/H] ^ —3 to —4). These results 
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are discussed in comparison with the previous observational studies along with the 
theoretical predictions from the available chemical evolution models. 



Key words: Galaxy: evolution — line: formation — stars: abundances — stars: atmospheres 
— stars: late-type 

1 INTRODUCTION 

The purpose of this paper is to clarify the behaviors of the abundances of sodium in disk/halo 
stars in our Galaxy. 

Recently, wc carried out an extensive non-LTE abundance analysis of potassium based on 
the K I 7699 resonance line and found a tendency of mildly supcrsolar [K/Fc] ratios for metal- 
poor stars just similar to the case of a-capturc elements (Takeda ct al. 2002). As a continuation 
of that study, wc pay attention this time to another similar alkali element, sodium (Na). 

The main site of Na synthesis is considered to be the hydrostatic carbon burning inside 
massive stars (e.g., Woosley, Weaver 1995), though other manufacturing processes (e.g., NeNa 
cycle in the H-burning shell, type I supernova, etc.) are possible and may have an appreciable 
affection on its chemical history (sec, e.g., Timmcs ct al. f995). As an oAA-Z (neutron- rich) 
clement, its production is sensitive to the neutron excess and therefore metal-dependent. Hence, 
clarifying the behavior of its abundance in terms of the cosmic mctallicity would provide us 
important information concerning the variation of the yield, the detailed mechanism of the 
synthesis, and the Galaxy evolution models (initial mass function, infall, etc ). 

Unlike the case of potassium, a reasonable number of adequate lines are available for sodium 
abundance determinations, which lie mostly in the visible through near-IR wavelength regions. 
Hence, many observational studies have so far been performed over these two decades (see, e.g., 
the references quoted in the recent theoretical papers such as Timmcs et al. 1995, Samland 
1998, and Goswami, Prantzos 2000). Nevertheless, the behavior of [Na/Fe] is currently still 
controversial and has not yet been settled. 

It has generally been considered until mid-1990's that Na almost scales with Fe at all 
metallicities though with a considerable scatter aroimd [Na/Fe] ~ 0. [See, for example, subsub- 
section 3.5.3 of Timmcs et al. (1995) for a review on the situation at the time of early 1990's.] 
Actually, a nearly solar ratio was obtained even for extremely metal-poor stars ([Fe/H] down 
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to ~ —4) by the extensive study of McWilliam et al. (1995b) using the resonance D hnes (Na i 
5890/5896). It should be pointed out that most of these previous studies were done based on 
the assumption of LTE. 

In late 1990's, however, new studies have appeared that cast doubt to this classical picture. 
Baumiiller et al. (1998) investigated the non-LTE formation of various Na i lines and found 
considerable (negative) non-LTE corrections reaching —0.5 dex for the D lines, by which a sys- 
tematic tendency of decreasing [Na/Fe] ratios is obtained toward a metal-poor regime ([Na/Fe] 

0.5 at [Fe/H] ~ -3). Similarly, the (LTE) analysis of Stephens (1999) for halo dwarf 

stars based on high-quality spectra of Na i 5683/5688 subordinate lines suggested that [Na/Fe] 
progressively decreases (from [Fe/H] —1 to [Fe/H] —2) down to the value of ~ —0.5. 

Accordingly, it appears that a proper account of the non-LTE effect is mandatory for 
studying stellar sodium abundances, at least for the resonance Na i 5890/5896 D lines visible 
even in extremely metal-poor stars. Then, what about other subordinate lines, such as Na i 
6154/6161 lines (widely used for relatively metal-rich population I disk stars), 5683/5688 lines 
(often used for typical halo stars with [Fe/H] of —1 to —2), 8183/8195 lines (not so popular due 
to their unfavorable wavelength location but potentially useful abundance indicators because 
of the clear visibility even down to [Na/H] ~ —3)? Do they also suffer an appreciable non-LTE 
effect ? Any definite conclusion on the Galactic [Na/Fe] vs. [Fe/H] behavior should await until 
this problem is clarified, because different lines are used case by case depending on the stellar 
metallicity. 

Admittedly, several elaborate and important studies have recently been done concerning 
the non-LTE effect on Na abundance determinations for solar- type stars (e.g., Baumiiller et 
al. 1998; Gratton et al. 1999; Korotin, Mishenina 1999; Mashonkina et al. 2000). We would 
say, however, that their results are not presented in sufficient detail for the purpose of practical 
applications; i.e., they tend to be confined to only specific lines of interest or only representative 
atmospheric parameters, and it is not necessarily easy for the reader to assess how all those 
Na I lines of importance are affected by the non-LTE effect for stars showing a diversity of 
atmospheric parameters. 

In view of this necessity, we decided to (1) carry out non-LTE calculations on the neutral 
sodium for a wide range of atmospheric parameters with an intention of applying to those 4 
line pairs mentioned above, (2) construct useful grids of non-LTE corrections, and (3) perform 
extensive abundance analysis toward establishing the Na abundances of disk/halo stars with 
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a variety of metallicities based on our own new observations/measurements along with the 
published equivalent widths taken from various literature. 

The following sections of this paper are organized as follows. Our statistical-equilibrium 
calculations on neutral sodium and the resulting non-LTE corrections are explained in section 2, 
where we also discuss uncertain factors or adopted approximations involved in the abundance 
determination based on some theoretical test calculations. Our observational data obtained 
at two observatories are explained in section 3, followed by section 4 where we describe the 
procedures of the abundance analysis using these data along with the reanalysis of the extensive 
literature data. We discuss in section 5 the results obtained in section 4, especially in terms 
of the [Na/Fe] vs. [Fe/H] relation and its implications, while referring to the representative 
theoretical calculations of Galactic chemical evolution. Finally, the conclusion of this paper is 
presented in section 6. 

2 NON-LTE CALCULATIONS 
2.1 Computational Procedures 

The procedures of our statistical-equilibrium calculations for neutral sodium, based on a Na I 
atomic model comprising 92 terms and 178 radiative transitions, arc essentially the same as 
those described in Takeda and Takada-Hidai (1994) and Takeda (1995), which should be con- 
sulted for details. One difFcrcncc is the choice of the background model atmospheres; i.e., instead 
of Kurucz's (1979) ATLAS6 models adopted in those old works, the photoionizing radiation 
field was computed based on Kurucz's (1993a) ATLAS9 model atmospheres while incorporating 
the line opacity with the help of Kurucz's (1993b) opacity distribution functions. It should also 
be mentioned that a reduction factor of 0.1 (logarithmic correction of /i = —1 according to 
the definition of Takeda 1995) was applied to the H i collision rates computed with classical 
approximate formula (Stecnbock, Holweger 1984; Takeda 1991) according to the empirical de- 
termination of Takeda (1995), though test calculations by varying this correction factor from 1 
to 10"^ were also performed (cf. subsection 2.3 below). 

Since we planned to make our calculations applicable to stars from near-solar metallicity 
(population I) down to very low metallicity (extreme population II) at late-F through early-K 
spectral types in various evolutionary stages (i.e., dwarfs, subgiants, giants, and supcrgiants), 
we carried out non-LTE calculations on an extensive grid of 125 (5x5x5) model atmospheres 
resulting from combinations of five Tcff values (4500, 5000, 5500, 6000, 6500 K), five log g values 
(1.0, 2.0, 3.0, 4.0, 5.0), and five metallicities (represented by [Fe/H]) (0.0, -1.0, -2.0, -3.0, 
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—4.0). As for the stellar model atmospheres, we adopted Kurucz's (1993a) ATLAS9 models 
corresponding to a microturbulent velocity (^) of 2 km s~^. 

Regarding the sodium abundance used as an input value in non-LTE calculations, we 
assumed logeJJ^"* = 6.33 + [Fe/H], where the solar sodium abundance of 6.33 was adopted 
from Anders and Grevesse (1989) (which is used also in the ATLAS9 models). Namely, a 
metallicity-scaled sodium abundance was assigned to metal-poor models. The microturbulent 
velocity (appearing in the line-opacity calculations along with the abundance) was assumed to 
be 2 km s~^, to make it consistent with the model atmosphere. 

2.2 Characteristics of the Non-LTE Effect 

In figure 1 are shown the Si,(t)/ B{t) (the ratio of the line source function to the Planck 
function, and nearly equal to ~ bu/b\, where bi and b^ are the non-LTE departure coefficients 
for the lower and upper levels, respectively) and Iq^^^ (t) / Iq^^ (t) (the NLTE-to-LTE line- 
center opacity ratio, and nearly equal to ~ b\) for each of the multiplet 1 (5890/5896) and 
multiplet 4 (8183/8195) transitions (non-LTE effects are especially important for these two 
doublets; cf. subsection 5.1) for a representative set of model atmospheres. We can read the 
following characteristics from this figure, which are mostly the same as those obtained for the 
case of K i 7699 (Takeda et al. 2002): 

— In almost all cases, the inequality relations oi Sl/B < 1 (dilution of line source function) 
and > 1 (enhanced line-opacity) hold in the important line-forming region for both 

cases of multiplets 1 and 4, which means that the non-LTE effect almost always acts in the 
direction of strengthening the 5890/5896 and 8183/8195 lines.* Actually, more important is 
the former S'L-dilution effect which appreciably deepens/darkens the core of saturated lines in 
such a way that would never bo accomplished under the assumption of LTE [where the residual 
intensity can not become lower than B{t ^ ())/ B{t ^ 1)]; hence, this raises the flat part of the 
curve of growth, and if this increase in the saturated-line strength is to be accounted for within 
the framework of LTE, a large abundance variation would have to be invoked. As a result, 
the extent of the non-LTE effect significantly depends on the line-strength in the sense that it 
becomes most conspicuous for the lines in the flat part. (See, e.g., Stiirenburg, Holweger 1990, 
especially figure 15 of subsection 4.3 therein.) This situation is demonstrated in figure 2, where 
the non-LTE corrections (see subsection 2.3 below) for the 5896 and 8195 lines are shown as 

* An exception is the case of lowest T^g / log g (e.g., Tgff =4500 K, log g = 1.0, [Fe/H] = 0.0), where the line can be 
marginally weakened (i.e., positive abundance correction) by the enhanced Sl over B in the line-forming region 
(cf. the footnote in subsection 5.3). 
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functions of the equivalent width for representative model atmospheres. 

— There is a tendency that the non-LTE effect is enhanced with a lowering of the gravity, as 
naturally expected. 

— The departure from LTE appears to be larger for higher Tes in the high-metallicity (Ix) 
case, while this trend becomes ambiguous, or even inverse, in the low-metallicity case. 

— Toward a lower metallicity, the extent of the non-LTE departure tends to decrease, but the 
departure appears to penetrate deeper in the atmosphere, which makes the situation rather 
complex. 

— For a very strong damping-dominated case (i.e., lowest T^a and highest metallicity), the 
departure from LTE shifts toward the upper atmosphere and the non-LTE effect becomes 
comparatively insignificant. 

2.3 Grid of Non-LTE Corrections 

Based on the results of these calculations, we computed extensive grids of the theoretical 
equivalent-widths and the corresponding non-LTE corrections for the considered eight lines 
(Na I 5683, 5688, 5890, 5896, 6154, 6161, 8183, and 8195) for each of the model atmospheres 
as follows. 

For an assigned sodium abundance {A'^) and microturbulcnce (^'*), we first calculated the 
non-LTE equivalent width (W^^"^^) of the line by using the computed non-LTE departure 
coefficients (6) for each model atmosphere. Next, the LTE {A^) and NLTE {A^) abundances 
were computed from this VK^'""'"^ while regarding it as if being a given observed equivalent 
width. We can then obtain the non-LTE abundance correction. A, which is defined in terms of 
these two abundances as A = A^ — A^ . 

Strictly speaking, the departure coefficients [&(t)] for a model atmosphere arc applicable 
only to the case of the sodium abimdance (log ej^'^'"'') and the microturbulcnce (2 km s^^) 
adopted in the non-LTE calculations (cf. subsection 2.1). Nevertheless, considering the fact 
that the departure coefficients (i.e., ratios of NLTE to LTE number populations) arc (unlike 
the population itself) not much sensitive to small changes in atmospheric parameters, we also 
applied such computed b values to evaluating A for slightly different and from those 
fiducial values assumed in the statistical equilibrium calculations. Hence, we evaluated A for 
three A'^ values (logej^'^'"' and ±0.3 dex perturbation) as well as three values (2 km s^^ and 
±1 km s^^ perturbation) for a model atmosphere using the same departure coefficients. 

We used the WIDTH9 program (Kurucz 1993a), which had been modified to incorporate 
the non-LTE departure in the line source function as well as in the line opacity, for calculating 



Non-LTE Analysis for the Sodium Abundances of Metal-Poor Stars in the Galactic Disk and Halo 



7 




Figure 1 Ratio of the line source function (S'l) to the local Planck func- 
tion (B) and the NLTE-to-LTE line-center opacity ratio as functions of the 
standard continuum optical depth at 5000 A computed for models of Tgg = 
4500 K, 5500 K, and 6500 K. Note that Si^/B and Z^l^e/^lte g^^.^ ^^^^^ -^^ 
the same line-type, but both are easily discernible since the former (diluted) 
is generally lower than the latter (overpopulated) . The solid lines show the 
results for the 3 ^S-3p transition of multiplet 1 (corresponding to Na i 
5890/5896 resonance lines), while those for the 3p ^P°-3d transition of 
multiplet 4 (corresponding to Na I 8183/8195 lines) are depicted by dashed 
lines. In each case, the results for two different gravity atmospheres are 
given: The thick lines are for logg — 4 and the thin lines are for log(? = 2, 
respectively. Note also that the curves are vertically offset by an amount of 
0.5 dex relative to those of the adjacent metallicity ones. 
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Figure 2 Theoretical non-LTE corrections as functions of the NLTE equiv- 
alent width for representative models, based on the electronic tables El (see 
subsection 2.3 for details). Different symbols discern the effective tempera- 
ture (triangle, inverse triangle, circle, square, and diamond correspond to Tes 
= 4500, 5000, 5500, 6000, and 6500 K, respectively) and the surface gravity 
(open and closed symbols are for log g = 2 and 4, respectively) of the adopted 
model atmospheres. The results derived for five different model-metallicities 
([Fe/H] = 0, —1, —2, —3, and —4), while using the metallicity-scaled sodium 
abundance along with perturbations of ±0.3 dex, are shown here (though 
the same symbols arc used for different metallicities/abundances). (a) Na I 
5895.92. (b) Na i 8194.82. 
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the equivalent width for a given abundance, or inversely evaluating the abundance for an as- 
signed equivalent width . The adopted line data {gf values, radiation damping constants, etc.) 
are given in table 1. 



Table 1 Adopted atomic data of the Na i lines. 



Mult. 


Line 


A 


X 


logs/ 


Tr 


1 


A5890 


5889.95 


0.00 


+0.12 


0.62 


1 


A5896 


5895.92 


0.00 


-0.18 


0.62 


4 


A8183 


8183.26 


2.10 


+0.22 


1.13 


4 


A8195 


8194.82 


2.10 


+0.52 


1.13 


5 


A6154 


6154.23 


2.10 


-1.56 


0.75 


5 


A6161 


6160.75 


2.10 


-1.26 


0.75 


6 


A5683 


5682.63 


2.10 


-0.67 


0.81 


6 


A5688 


5688.21 


2.10 


-0.37 


0.81 



Notes. Each of the columns 1 6 give the multiplet number, the line abbreviation used in this 
paper, the line wavelength (in A), the lower excitation potential (in eV) , the logarithm of the 
gf value, and the radiation damping constant in unit of 10^ s~^. The data presented here are 
based on table 1 of Takeda and Takada-Hidai (1994), who consulted the compilation of Wiese 
et al. (1969). Each line was treated as if it is a single line, while neglecting any hyperfine 
structure (cf. subsection 2.3). See also subsection 2.3 for the treatment of the quadratic Stark 
effect damping and of the van der Waals effect damping. 



One of the controversial factors in abundance determinations for late-type stars is the 
choice of the van der Waals effect damping constant. Regarding this parameter, we assumed 
the classical Unsold's (1955) formula imchanged, which means the adoption of AlogCg = 0.0 
(Cg is connected to the damping width as logFe = logrg'''^*'"^''' + 0.4AlogC6). While this choice 
is based on our previous empirical investigations (cf. appendix A of Takeda, Takada-Hidai 1994; 
subsection 4.2 of Takeda 1995), we actually confirmed that this choice is reasonable from the 
analysis of solar Na i lines (cf. subsection 4.2). Anyhow, the precise value of this correction 
is not very essential unless very strong lines are used, as can be seen from figure 3 where the 
abundance changes expected by using AlogCg = 0.5 are graphically displayed. Namely, the 
maximum change is ^ 0.4|A logCej at most occurring only in the case of very strong damping- 
dominated lines and the low-temperature condition where Eg dominates E (= Er + E4 + Eg). 
Regarding the quadratic Stark effect damping (E4, which is comparatively insignificant in late- 
type stars), we followed the Feytremann's formula (Kurucz 1979, p. 8). 

Note also that we neglected the hyperfine structure (hfs; see, e.g., table 4 of McWilliam et 
al. 1995b for the 5890/5896 hues, and table 1 of Takeda, Takada-Hidai 1994 for the 8195 line) 
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Figure 3 Differences of the resulting non-LTE abundances, when 
AlogCe = +0.5 (logarithmic correction applied to classical value of the 
van der Waals effect damping parameter) was used instead of the fiducial 
value of AlogCg = 0.0, as functions of the non-LTE equivalent widths for 
the representative cases shown in table 2. See the caption of figure 2 for the 
meanings of the symbols, (a) Na i 5895.92. (b) Na i 8194.82. 
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while treating each line as being purely single, because of a technical reason.''' However, this 
does not make any serious affection (several hundredths dex in most cases of actual importance) 
as illustrated in figure 4. 

As a demonstrative example of non-LTE corrections, we give the ^ = 1 km s~^ results for the 
Na I 5896 and 8195 lines computed for representative parameters [chosen to be nearly compatible 
with Baumiiller et al.'s (1998) calculations] in table 2, where we also present the cases oih = 
0, —2, —3 in addition to the fiducial h = —1 case {h is the logarithmic H i collision correction) 
for a comparison. Comparing this table 2 with Baumiiller et al.'s (1998) table 2 (where they 
adopted lO'' = 0.05 or /i = —1.3), we can see that our |A| values are quantitatively somewhat 
larger than theirs, though the qualitative tendency is in good agreement. The complete results 
(for all combinations of Tefj, logg, ^ values for each of the 8 lines, though only for the fiducial 
case oih= —1) are available only in the electronic form. These tables (named "tables El") are 
temporarily at the anonymous ftp site of 

\protect\vrule widthOpt\protect\href {f tp: //www. ioa. s .u-tokyo.ac . jp/Users/takeda/sodiuiii_paper/}-{f tp: //www. ioa. 
(IP address: 133.11.160.242) 

but will be replaced to CDS when the paper is to be published. 
3 OBSERVATIONAL DATA 

Our observational data consist of those obtained at two observatories, Beijing Astronomical 
Observatory (BAO; 2.2 m reflector + coude echelle spectrograph at Xinglong station) in P. R. 
China and Okayama Astrophysical Observatory (OAO; 1.9 m reflector + coude High-Dispersion 
Echelle Spectrograph named "HIDES" ) in Japan. 

3.1 BAO Data 

The BAO spectra {R ^ 40000 and S/N ~ 200-400) used here are such those originally collected 
for the purpose of Chen et al.'s (2000) comprehensive analyses on Galactic disk stars. See 
subsection 2.2 therein for detailed information on the observations and the data quality. Chen 
et al. (2000) determined the abundance of sodium based on the comparatively weak Na i 
6154/6161 lines, though the used equivalent- widths data were not published. For the present 

tWe decided to invoke the original WIDTH9 algorithm (applicable to a symmetric single line) for an equivalent- 
width calculation instead of using the spectral-synthesis technique (which requires a sufhciently fine division of 
a line profile unsuited / unnecessary for very strong lines) , since we wanted to treat cases of considerably different 
line strengths in the same consistent way (i.e., for constructing the tables of non-LTE corrections for different 
lines over wide parameter ranges). 
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Figure 4 Differences of the resulting non-LTE abundances, when the hy- 
perfine structure was approximately included instead of the purely single-line 
treatment (as basically adopted in this study; cf. table 1), as functions of 
the non-LTE equivalent widths for the representative cases shown in table 
2. See the caption of figure 2 for the meanings of the symbols, (a) Na i 
A5896 [using the two-component approximation of (//(5895.93) = 0.269 and 
p/(5895.95) = 0.376; cf. McWilham et al. 1995b]. (b) Na i A8195 [using the 
two-component approximation of 5/(8194.79) = 0.332 and 5-/(8194.82) = 
2.999; cf. Wiese et al. 1969]. 
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Table 2 Dependence of the non-LTE effect on the choice of the correction 
factor for the H l collision. 





log 3 


[Fc/HJ 


c 


A" 


A 


Wi, 


W-i 


W-2 


W-3 


Ao 


A-i 


A-2 


A-3 


5500 


4.0 


0.0 


1.0 


6.33 


5895.92 


602.6 


616.6 


616.6 


631.0 


-0.04 


-0.06 


-0.06 


-0.06 


5500 


4.0 


-1.0 


1.0 


5.33 


5895.92 


309.0 


316.2 


331.1 


331.1 


-0.10 


-0.15 


-0.17 


-0.18 


5500 


4.0 


-2.0 


1.0 


4.33 


5895.92 


158.5 


173.8 


177.8 


177.8 


-0.26 


-0.37 


-0.40 


-0.41 


5500 


4.0 


-3.0 


1.0 


3.33 


5895.92 


85.1 


89.1 


91.2 


89.1 


-0.24 


-0.33 


-0.35 


-0.33 


5500 


4.0 


-4.0 


1.0 


2.33 


5895.92 


20.0 


20.9 


21.4 


20.9 


-0.03 


-0.06 


-0.06 


-0.05 


6500 


4.0 


0.0 


1.0 


6.33 


5895.92 


263.0 


269.1 


275.4 


275.4 


-0.18 


-0.21 


-0.21 


-0.22 


6500 


4.0 


-1.0 


1.0 


5.33 


5895.92 


169.8 


177.8 


182.0 


182.0 


-0.47 


-0.55 


-0.56 


-0.56 


6500 


4.0 


-2.0 


1.0 


4.33 


5895.92 


109.7 


114.8 


114.8 


114.8 


-0.61 


-0.71 


-0.72 


-0.72 


6500 


4.0 


-3.0 


1.0 


3.33 


5895.92 


38.0 


39.8 


39.8 


38.9 


-0.13 


-0.16 


-0.16 


-0.16 


6500 


4.0 


-4.0 


1.0 


2.33 


5895.92 


5.1 


5.4 


5.4 


5.4 


-0.06 


-0.08 


-0.08 


-0.08 



5500 


4.0 


0.0 


1.0 


6.33 


8194.82 


295.1 


309.0 


316.2 


316.2 


-0.15 


-0.20 


-0.21 


-0.22 


5500 


4.0 


-1.0 


1.0 


5.33 


8194.82 


154.9 


173.8 


182.0 


186.2 


-0.20 


-0.35 


-0.42 


-0.44 


5500 


4.0 


-2.0 


1.0 


4.33 


8194.82 


57.5 


66.1 


69.2 


70.8 


-0.08 


-0.20 


-0.25 


-0.26 


5500 


4.0 


-3.0 


1.0 


3.33 


8194.82 


10.0 


11.8 


12.3 


12.3 


-0.03 


-0.10 


-0.13 


-0.13 


5500 


4.0 


-4.0 


1.0 


2.33 


8194.82 


1.1 


1.3 


1.4 


1.4 


-0.03 


-0.10 


-0.12 


-0.12 


6500 


4.0 


0.0 


1.0 


6.33 


8194.82 


199.5 


213.8 


213.8 


218.8 


-0.34 


-0.43 


-0.47 


-0.47 


6500 


4.0 


-1.0 


1.0 


5.33 


8194.82 


107.2 


117.5 


123.0 


123.0 


-0.38 


-0.55 


-0.60 


-0.61 


6500 


4.0 


-2.0 


1.0 


4.33 


8194.82 


34.7 


38.0 


38.9 


38.9 


-0.12 


-0.19 


-0.21 


-0.21 


6500 


4.0 


-3.0 


1.0 


3.33 


8194.82 


4.9 


5.5 


5.6 


5.6 


-0.09 


-0.13 


-0.14 


-0.14 


6500 


4.0 


-4.0 


1.0 


2.33 


8194.82 


0.5 


0.6 


0.6 


0.6 


-0.09 


-0.13 


-0.14 


-0.14 



Notes. Columns 1-6 are self-explanatory (the units of Tgg, g, and ^ are K, cm s ^, and km s ^, 

respectively). The suffixes (0, —1, —2. —3. and —4) appended to W (the non-LTE equivalent width in 
mA calculated for the atmospheric parameters and the assigned sodium abundance given in columns 
1-5) and A (the non-LTE abundance correction) denote the corresponding values of h (the logarithm 
of the H I collision correction factor applied to the classical formula) . 
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study, additional measurements for five lines (Na i 5683, 5688, 5890, 5896, and 8195; the 8183 
line was discarded because of the contamination of telluric lines) were newly carried out by 
one of us (Y.-Q. Chen; cf. subsection 2.3 of Chen et al. 2000 for the measurement method) 
on the spectra of 22 F-G stars (+ Moon), which are our BAG program stars also adopted 
in our previous potassium analysis (Takeda et al. 2002)*. Although many lines (except for 
the 6154/6161 lines) are so strong in these disk stars and not suitable for sodium abundance 
determinations, we tried to derive the abundances from them in order to see whether the similar 
behavior of [Na/Fe] may be reproduced among the different line groups (cf. subsection 5.2). 
Yet, it should be kept in mind that difficulties are involved in measuring the equivalent widths 
of very strong damping-dominated lines such as the 5890/5896 D lines, with which errors are 
more or less involved. Nevertheless, comparing the solar (Moon) equivalent width measured 
by us with those taken from the literature (cf. table 3), we can see that our measurements 
are regarded as being reasonable (i.e., crucial systematic discrepancies are not observed). The 
finally obtained BAG equivalent- widths are given in tables 4 and 5. 

3.2 GAG Data 

The GAG data are based on the near-IR spectra of 17 mostly metal-poor disk/halo stars (dwarfs 
and giants being mixed), which have been collected by using HIDES during our observing runs 
in 2001-2002.§ The slit width of 250 yum (0."95) was set to yield a spectral resolution of 
R ~ 50000. By using the single 4Kx2K CCD (pixel size of 13.5 /^m x 13.5 /xm), a wavelength 
span of ~ 1100 A could be covered at one time. Most of our observations were made in the 
wavelength region of ~ 7700-8800 A. 

The data reduction was performed with the IRAF^ echelle package, following the stan- 
dard procedure for extracting one-dimensional spectra. We also removed the telluric lines using 
the telluric task of IRAF and the spectra of a rapid rotator. The S/N ratios of the resulting 
spectra, differing from star to star, are typically of the order of 200-300. The resulting GAG 
equivalent widths of the Na i 8183/8195 lines based on these spectra, which were measured by 
the Gaussian fitting or direct integration method using the splot task of IRAF, are presented 



tOne of our present 22 BAO stars, HD 167588, was not included in Takeda et al.'s (2002) analysis, because there 
was a problem in the data quality of the K i 7699 line. 

§ Actually, the primary motivation of these observations was to study the abundance of sulfur based on the near- 
IR S I lines. Nevertheless, we could apply those spectra to the present purpose, since they cover the wavelength 
region of the Na i 8183/8195 lines. 

^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Tabic 3 Compaiison of the solar equivalent widths (in mA) of Na I lines 
taken from various literature. 



Ref. 


A5683 


A5688 


A5890 


A5896 


A6154 


A6161 


A8183 


A8195 


Data source 


(1) 


100.4 


131.1 


837.2 


545.0 


38.6 


59.5 




306.1 


BAO Moon" 


(2) 


105.3 


128.6 














B76 solar flux 


(3) 


106.9 


127.9 














B76 solar flux 


(4) 


106.9 


127.9 






36.8 


58.6 






B76 solar flux 


(5) 




129.0 


1125.5 


736.5 






257.3 


326.6 


solar flux & disk center'' 


(6) 










35.0 


53.7 






solar disk center'^ 


(7) 


90.0 


119.1 






40.5 


58.2 






K84 solar flux 


(8) 


121 


144 


830 


640 


41 


63 


254 


328 


K84 solar flux 


(9) 






756 


597 






244 


315 


K84 solar flux'' 


(10) 


103.0 


128.0 






36.0 








From (12)? 


(11) 










38.6 


59.5 






K84 solar flux 


(12) 


103 


128 


765 


570 


36 


53 


239 


322 


solar disk center^ 



References: (1) This study; (2) Peterson and Carney (1979); (3) Peterson (1980); (4) Peterson 
(1981); (5) Gratton and Sneden (1987b); (6) Tomkin et al. (1985); (7) Prochaska et al. (2000); 
(8) Mashonkina et al. (2000): (9) Takeda (1995); (10) Carretta et al. (2000); (11) Sadakane et 
al. (2002); (12) Holweger (1971). 



Notes on the data sources: 

B76 and K84 means the solar flux spectrum atlas published by Beckers et al. (1976) and Kurucz 
et al. (1984), respectively. 

Moon spectrum observed at Beijing Astronomical Observatory. 

Based on two different atlases (see Gratton, Sneden 1987a). 
" Taken from the values published by Lambert and Luck (1978), which are based on two dif- 
ferent atlases. 

Inversely computed from the profile fitting solutions. 
^ Based on three different atlases. 
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Table 4 Analysis of Na i 8183/8195 lines based on the BAO and OAO data. 



Sla. 




Sp. lypc 


■J off 


log</ 


[I'r/'ll^ 


S 








11 ^lor, 




A.si,,r, 


lirl'. 












[BAO sample] 
















Sun 




G2 V 


5780 


4.44 


0.00 


1.0 








306.1 


6.32 


-0.19 


(1) 


HD 


010307 


G1.5 V 


5776 


4.13 


-0.05 


1.8 








265.5 


6.08 


-0.31 


(1) 


HD 


019373 


GO V 


5867 


4.01 


+0.03 


1.8 








274.8 


6.23 


-0.32 


(1) 


HD 


022484 


F9 IV-V 


5915 


4.03 


-0.13 


2.0 








242.5 


5.97 


-0.39 


(1) 


HD 


034411 


G1.5 IV-V 


5773 


4.02 


+0.01 


1.7 








275.8 


6.20 


-0.30 


(1) 


HD 


039587 


GO V 


5805 


4.29 


-0.18 


2.2 








273.9 


5.99 


-0.32 


(1) 


HD 


041640 


F5 


6004 


4.37 


-0.62 


2.0 








205.6 


5.65 


-0.39 


(1) 


HD 


049732 


F8 


6260 


4.15 


-0.70 


1.9 








178.3 


5.63 


-0.50 


(1) 


HD 


055575 


GO V 


5802 


4.36 


-0.36 


1.6 














(1) 


HD 


060319 


F8 


5867 


4.24 


-0.85 


1.6 








171.4 


5.41 


-0.39 


(1) 


HD 


062301 


F8 V 


5837 


4.23 


-0.67 


1.7 








190.1 


5.53 


-0.38 


(1) 


HD 


068146 


F7 V 


6227 


4.16 


-0.09 


2.1 








229.4 


6.01 


— 0.42 


(1) 


HD 


069897 


F6 V 


6243 


4.28 


-0.28 


2.0 








197.5 


5.77 


-0.42 


(1) 


HD 


076349 


GO 


6004 


4.21 


-0.49 


2.1 








210.9 


5.69 


-0.43 


(1) 


HD 


101676 


F6 V 


6102 


4.09 


-0.47 


2.0 








214.8 


5.82 


-0.47 


(1) 


HD 


106516 


F5 V 


6135 


4.34 


-0.71 


1.5 








168.6 


5.55 


-0.41 


(1) 


HD 


109303 


F8 


5905 


4.10 


-0.61 


1.7 








186.5 


5.56 


-0.43 


(1) 


HD 


118244 


F5 V 


6234 


4.13 


-0.55 


2.3 








197.7 


5.68 


-0.50 


(1) 


HD 


142373 


F8 Ve 


5920 


4.27 


-0.39 


1.5 














(1) 


HD 


142860 


F6 IV 


6227 


4.18 


-0.22 


2.2 








211.5 


5.84 


-0.44 


(1) 


HD 


167588 


F8 V 


5894 


4.13 


-0.33 


1.7 














(1) 


HD 


201891 


F8 V-IV 


5827 


4.43 


-1.04 


1.6 








152.6 


5.21 


-0.33 


(1) 


HD 


208906 


F8 V-IV 


5929 


4.39 


-0.73 


1.5 








173.9 


5.47 


-0.36 


(1) 












[OAO sample] 
















HD 


6833 


G9 III 


4450 


1.40 


-0.90 


1.6 


161.3 


4.92 


-0.41 


207.9 


5.10 


-0.46 




HD 


26297 


G5/G6 Ivw 


4500 


1.20 


-1.60 


1.7 


86.4 


4.24 


-0.21 


115.6 


4.24 


-0.32 


l^i 


HD 


73394 


G5 IIIw 


4500 


1.10 


-1.40 


1.5 


127.7 


4.71 


-0.37 


157.0 


4.73 


-0.46 


(3) 
(4)" 


HD 


76932 


F7/F8 IV/V 5900 


4.12 


-0.80 


1.3 


125.1 


5.40 


-0.38 


155.6 


5.37 


-0.42 


HD 


88609 


G5 IIIw 


4570 


0.75 


-2.70 


1.9 


10.8 


3.18 


-0.10 


20.1 


3.19 


-0.11 


(4) 
(4)" 


HD 


106516 


F5 V 


6200 


4.31 


-0.70 


1.1 


134.8 


5.66 


-0.41 


167.1 


5.65 


-0.42 


HD 


108317 


GO 


5300 


2.90 


-2.24 


1.0 


26.5 


3.95 


-0.16 


41.0 


3.91 


-0.18 


(5) 


HD 


122563 


F8 IV 


4590 


1.17 


-2.74 


2.3 


15.7 


3.35 


-0.11 


31.5 


3.41 


-0.12 


(6) 
(7)° 


HD 


140283 


sdF3 


5690 


3.69 


-2.42 


0.8 


6.0 


3.38 


-0.12 


9.6 


3.29 


-0.13 


HD 


165908 


F7 V 


5900 


4.09 


-0.60 


1.7 


143.2 


5.51 


-0.39 


179.9 


5.51 


-0.43 


(4) 


HD 


167588 


F8 V 


5890 


4.13 


-0.33 


1.7 


170.9 


5.76 


-0.38 


205.3 


5.72 


-0.39 


(1) 


HD 


187111 


GSwvar 


4260 


0.51 


-1.85 


1.8 


94.5 


4.18 


-0.19 


147.6 


4.41 


-0.34 


(6) 


HD 


189322 


G8 III 


4464 


2.00 


-1.57 


2.0 


222.1 


5.37 


-0.43 


260.0 


5.38 


-0.38 


(8)" 


HD 


216143 


G5 


4525 


1.00 


-2.10 


2.9 


50.0 


3.86 


-0.12 


70.7 


3.78 


-0.15 


(2) 


HD 


221170 


G2 IV 


4425 


1.00 


-2.15 


1.5 


53.0 


3.91 


-0.14 


79.4 


3.93 


-0.20 


(5) 


BD 


+37° 1458 GO 


5200 


3.00 


-2.00 


1.3 


19.8 


3.73 


-0.15 


42.7 


3.86 


-0.18 


(2) 


Procyon 


F5 IV-V 


6510 


3.96 


-0.05 


2.2 


188.7 


6.14 


-0.51 


220.3 


6.08 


-0.53 


(9) 



Notes. Columns 1-6 are self-explanatory as in tabic 2. W is the observed equivalent width (in mA), is the 

logarithmic non-LTE abundance (in the usual normalization of H — 12), and A is the non- LTE correction (= A" - A""), 
respectively. Cohimn 13 gives the keys to the references of the atmospheric parameters (cf. subsection 4.1): (1) Chen 
ct al. (2000);(2) Fulbrighl (2000): (3) Luck and Bond (1985);(4) Takada-Hidai ct al. (2002): (5) Pilachowski et al. 
(1996); (6) Gratton and Sneden (1994); (7) Nissen et al. (2002); (8) Alonso et al. (1999); (9) AUende Prieto et al. 
(2002). 

Only the ^ value was taken from (2). 
^ The ^ value was assumed. 
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Table 5 Analysis of Na i 5683/5688, 5890/5896, and 6154/6160 lines based 
on the BAO data. 





W5683 ^S6R3 


A 5683 


W5688 ^5688 ^^5688 


^5890 ^5890 


^5890 


-rrz 

VV5896 




T 

-^5896 


V> 6154 


^6154 ^ 


\^ 

-^6 


154 


V> 6161 


4'' 
f? 


61 -^6 


61 
























[BAO sample! 


























Sun 


100.4 


6 


28 


-0 


11 


131.1 


6 


33 


-0 


13 


837.2 


6.41 


-0.05 


545.0 


6.27 


— 0.07 


38.6 


6.31 


— 


04 


59.5 


6 


32 


-0 


06 


xjT~v ninQn'7 












134.3 


6 


26 


-0 


16 


779.2 


6.43 


-0.07 


506.5 


6.30 


—0.09 


39.4 


6.29 


— 


05 


60.6 


6 


28 


-0 


07 


HD 019373 












139.6 


6 


39 


-0 


19 














45.8 


6.44 


-0 


05 


65.3 


6 


40 


-0 


07 


HD 022484 












116.9 


6 


12 


-0 


17 


458.5 


6.02 


-0.12 


402.6 


6.16 


-0.13 


32.1 


6.23 


-0 


06 


47.1 


6 


17 


-0 


07 


HD 034411 












133.6 


6 


30 


-0 


17 


795.3 


6.52 


-0.07 


505.6 


6.37 


-0.08 


39.7 


6.30 


-0 


05 


65.6 


6 


36 


-0 


07 


HD 039587 












118.1 


6 


01 


-0 


13 


625.6 


6.11 


-0.09 


459.5 


6.06 


-0.12 


24.0 


6.01 


-0 


05 


48.4 


6 


11 


-0 


06 


HD 041640 












77.2 


5 


69 


-0 


12 


413.1 


5.69 


-0.16 


316.3 


5.64 


-0.21 


9.1 


5.63 


-0 


07 


21.8 


5 


75 


-0 


07 


HD 049732 












73.7 


5 


78 


-0 


15 


343.8 


5.80 


-0.24 


260.8 


5.65 


-0.34 


7.7 


5.65 


-0 


08 


17.3 


5 


74 


-0 


09 


HD 055575 












93.2 


5 


82 


-0 


12 








412.6 


5.90 


-0.11 


21.6 


5.96 


-0 


05 


36.2 


5 


95 


-0 


06 


HD 060319 












59.7 


5 


45 


-0 


11 


354.5 


5.42 


-0.18 


270.9 


5.36 


-0.23 


6.3 


5.40 


-0 


07 


17.8 


5 


60 


-0 


08 


HD 062301 












71.8 


5 


58 


-0 


12 


407.8 


5.60 


-0.14 


285.8 


5.43 


-0.21 


10.2 


5.61 


-0 


07 


19.9 


5 


64 


-0 


07 


HD 068146 












115.9 


6 


24 


-0 


17 


436.3 


6.19 


-0.15 


352.0 


6.20 


-0.18 


24.6 


6.23 


-0 


06 


43.3 


6 


25 


-0 


07 


HD 069897 












94.1 


6 


01 


-0 


14 


374.0 


5.92 


-0.18 


306.3 


5.92 


-0.23 


14.5 


5.96 


-0 


06 


29.1 


6 


02 


-0 


07 


HD 076349 












87.1 


5 


79 


~0 


14 


401.0 


5.74 


-0.17 


318.0 


5.71 


-0.22 


14.7 


5.85 


-0 


07 


29.8 


5 


92 


-0 


08 


HD 101676 












93.0 


5 


93 


-0 


16 








294.5 


5.76 


-0.27 


15.3 


5.92 


-0 


07 


31.2 


5 


99 


-0 


08 


HD 106516 












69.9 


5 


71 


-0 


13 


313.5 


5.52 


-0.21 


260.6 


5.56 


-0.26 


9.0 


5.68 


-0 


07 


22.1 


5 


82 


-0 


08 


HD 109303 












78.6 


5 


70 


-0 


14 


382.4 


5.66 


-0.17 


279.4 


5.53 


-0.23 


8.1 


5.53 


-0 


07 


32.2 


5 


93 


-0 


08 


HD 118244 












91.9 


5 


93 


-0 


16 


353.7 


5.78 


-0.25 


281.9 


5.68 


-0.34 


9.7 


5.75 


-0 


08 


29.5 


6 


00 


-0 


09 


HD 142373 


68.1 


5 


89 


-0 


10 












407.0 


5.75 


-0.12 


401.1 


6.03 


-0.13 


12.6 


5.75 


-0 


06 


26.5 


5 


84 


-0 


07 


HD 142860 












106.3 


6 


11 


-0 


16 


449.6 


6.18 


-0.16 


317.6 


5.98 


-0.24 


17.6 


6.05 


-0 


06 


44.0 


6 


25 


-0 


08 


HD 167588 


75.1 


5 


96 


-0 


11 


97.4 


5 


92 


-0 


15 


423.4 


5.84 


-0.13 


368.4 


5.96 


-0.14 


19.1 


5.95 


-0 


06 


37.6 


6 


02 


-0 


07 


HD 201891 












43.3 


5 


20 


-0 


09 


328.1 


5.14 


-0.19 


209.3 


4.82 


-0.28 


7.4 


5.46 


-0 


07 


9.6 


5 


28 


-0 


07 


HD 208906 


52.1 


5 


68 


-0 


09 
























9.2 


5.61 


-0 


07 


21.3 


5 


72 


-0 


07 



See the notes to table 4 for the meanings of the W, A , and A. 



in tabic 4. Note that, although two stars (HD 106516 and HD 167588) are common to BAO 
and OAO samples, they were treated independently. 



4 ABUNDANCE ANALYSES 
4.1 Atmospheric Parameters 

First we need to establish the four atmospheric parameters [Tcs (effective temperature), logg 
(surface gravity), [Fe/H] (mctaUicity; represented by the Fc abundance relative to the Sun), 
and ^ (microturbulcnt velocity dispersion)], which arc necessary for constructing the model 
atmosphere and deriving the abundance from an equivalent width. Regarding the 22 BAO 
stars, the parameter data determined by Chen ct al. (2000) were used unchanged, while the 
reasonably known solar parameters of (5780 K, 4.44, 0.0, 1.0 km s~^) were applied to the 
analysis of the BAO Moon data. Meanwhile, the atmospheric parameters of 17 OAO stars were 
taken from the published values found in various literature (cf. column 13 of table 4); when two 
or more references were available for the same star, we selected an appropriate value according 
to our personal judgment. The finally adopted values are given in table 4. 

As for the model atmospheres, Kurucz's (1993a) grid of ATLAS9 models was used as in 
the case of non-LTE calculations, based on which the model of each star was obtained by a 
three-dimensional interpolation with respect to T^g, logg, and [Fe/H]. Similarly, the depth- 
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dependent departure coefficients (6) of Na i levels computed for the grid of models (cf. section 
2) were interpolated in terms of Tes, loggf, and [Fe/H], in order to evaluate the Si,{t)/B{t) and 
;NLTE(^y;LTE(^-) ratios for each star. 

4.2 Abundance Determination 

The procedures for determining the sodium abundances from the equivalent widths of Na I lines 
are the same as already explained in subsection 2.3 (i.e., modified WIDTH9 program, line data 
given in table 1, AlogCg = 0, neglecting the hfs effect). 

The resulting non- LTE abundance (A^) and the non-LTE abundance correction {— A — 
A^) derived for each line and for each star are presented in tables 4 (8183/8195 lines) and 5 
(5683/5688, 5890/5896, and 6154/6161 lines). As can be seen from these two tables, the solar 
non-LTE abundances derived from the BAO Moon equivalent widths of Na I 8195 (6.32), 5683 
(6.28), 5688 (6.33), 5890 (6.41), 6154 (6.31), and 6160 (6.32) are in remarkable agreement with 
each other; moreover, their average of 6.32((t = 0.04) is essentially the same as the standard 
solar sodium abundance of 6.33 (Anders, Grevesse 1989). Accordingly, we confirmed the in- 
ternal consistency of using AlogCg = 0.0, as far as our present analysis using ATLAS9 model 
atmospheres is concerned. [Actually, empirical determination of this correction significantly 
depends on the choice of the atmospheric models, as demonstrated in subsection 4.2 of Takeda 
(1995).] 

The abundance changes caused by uncertainties in the atmospheric parameters are quan- 
titatively similar to the case of potassium described in section 5 of Takeda et al. (2002), since 
K and Na have quite similar atomic structures to each other (alkali atom with one valence 
electron; ionization potential of 4- 5 eV) and most atoms are in the once-ionized stage. Hence, 
even for changes of ATcff ^ ±200 K or Alogg ^ ±0.3 dex (we expect that internal errors in 
the parameters of om program stars, especially for BAO sample stars, are smaller than these), 
the extents of the resulting abundance variations are ^0.2 dex, as can be seen from table 1 
of Takeda et al. (2002). The effect of changing the microturbulcnce, which influences only 
the lines at the flat part of the curve of growth {W\ ~ 100-300 mA; cf. figure 4) can also be 
assessed by inspecting that table. 

The [Na/Fe] ratios (Na-to-Fe logarithmic abundance ratio relative to the Sun) can be ob- 
tained based on the results in tables 4 and 5 as yl^(star) — 6.33 — [Fe/H], where we assumed 
Anders and Grevesse's (1989) value of 6.33 as the reference solar Na abundance. The resulting 
[Na/Fe] vs. [Fe/H] relations derived from each of the four multiplets, along with the corre- 
sponding non-LTE corrections, are displayed in figure 5. 
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Figure 5 The [Na/Fe] values resulting from the non-LTE analysis of our 

BAO and OAO data (larger symbols; scale in the left-hand axis) and the 
corresponding non-LTE corrections (smaller symbols; scale in the right-hand 
axis), plotted as functions of [Fe/H]; those arc based on the data presented 
in tables 4 and 5. (a) A5683/A5688 lines of multiplet 6, plotted in half- 
filled squares (BAO). (b) A5890/A5896 lines of multiplet 1, plotted in filled 
squares (BAO). (c) A6154/A6161 lines of multiplet 5, plotted in open squares 
(BAO). The results of our reanalysis of Edvardsson et al.'s (1993) data (also 
depicted in figure 7a) and the results of Feltzing and Gustafsson (1998) (the 
LTE [Na/Fe] values obtained by themselves; their equivalent-width data are 
not published) are also shown with small open circles and dots, respectively, 
for a comparison, (d) A8183/A8195 lines of multiplet 4, plotted in Greek 
crosses (BAO) and St. Andrew's crosses (OAO), respectively. 
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Note that when two Unes of the same multiplet are available, we averaged the abundances 
([Na/Fe]) as well as the non-LTE corrections for both lines and such derived mean results for 
the multiplet are shown. 

4.3 Literature Data Analysis 

We also tried to reanalyze the published equivalent-width data (for the Na I 5683/5688, 5890/5896, 
6154/6161, and 8183/8195 lines) of late-type stars in the Galactic disk as well as in the halo 
covering a wide range of mctallicities, while taking account of the non-LTE effect. In searching 
for the papers including the observational data we need, the extensive references quoted by 
Timmcs ct al. (1995), Samland (1998), and Goswami and Prantzos (2000) were quite helpful. 
Although our literature survey is not complete, we consider that we have picked up most of the 
important works done after 1980's. 

These data were analyzed just in the same way as we did for owe BAO/OAO data (cf. 
subsections 4.1 and 4.2). The atmospheric parameters were taken from the same paper as that 
presenting the equivalent-width data as far as possible." By inspecting those original papers, we 
foimd that abundance analyses done before 1990 tend to be based on rather coarsely determined 
parameters (e.g., using rounded Toff or log g values; assuming the same ^ values for all program 
stars, etc.) compared to the more recent studies carried out this decade. Accordingly, we 
decided to present the results of the reanalyzes on the data before and after 1990 separately. 
Figure 6 shows the resulting [Na/Fe] vs. [Fe/H] relations and the corresponding non-LTE 
corrections derived from the older data before 1990, while those obtained from the analyses of 
the recent data after 1990 arc displayed in figure 7.** 

The same results as those in figure 7a (the rcanalysis results based on the new data after 
1990) are also depicted in figure 8a, where the [Na/Fe] data arc plotted with a more compressed 
vertical scale (for the purpose of clarifying the global tendency) and low- and high-gravity stars 
are distinguished by different symbols. In addition, figure 8b shows the Tcff vs. logg plots 
relevant to these data; it can be seen from this figure that most of the high-gravity stars are 
comparatively higher Tcff stars of mid-F through late-G type, while low-gravity stars mainly 
consist of late-G through early-K giants of lower Teff and show a positive correlation between T^g 

II The exceptional cases are as follows: Gratton and Sneden's (1987b) EW data were analyzed with the parameters 
taken from Gratton and Snedcn (1987a). Peterson's (1989) EW data were analyzed with the parameters of 
Peterson et al. (1990). Zhao and Magain's (1990b) EW data were analyzed with (T^fi, logg, and ^) values taken 
from Magain (1989) and [Fe/H] values taken from Zhao and Magain(1990a). McWilliam et al.'s (1995a) EW 
data were analyzed with the parameters taken from McWilliam et al. (1995b). 

**Note that the BAG results are shown only for the 6156/6161 lines in figure 7, because those for the other strong 
lines are comparatively less reliable. 
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and log 5 (reflecting the evolutionary sequence). We may state from an inspection of figure 8a 
that any meaningful systematic log ^-dependent difference does not exist in the [Na/Fe] values; 
however, the biased sample (i.e., very low-metallicity stars tend to be low-gravity giants, and 
vice versa) prevents us from making any definite argument regarding this point. 

The details of these analyses (the data of the used equivalent widths and the adopted 
parameter values, the resulting non-LTE abundances or [Na/Fe] values with the non-LTE cor- 
rections, given for each line/multiplet and for each star) are available only electronically as 
"tables E2"(from the ftp site mentioned at the end of subsection 2.3). 

5 DISCUSSION 
5.1 Non-LTE Effect 

Figures 5, 6, and 7 clearly show how the non-LTE effect acts on different Na i lines in sodium 
abundance determinations for stars of various metallicities. 

Roughly speaking, a simple principle appears to hold: The non-LTE effect (negative correc- 
tions to the corresponding LTE abundances) tends to be most important for rather strong satu- 
rated lines (with equivalent widths of 100 300 mA) on the flat part of the curve of growth, and 
it becomes progressively insigniflcant towards decreasing strengths (i.e., imsaturated weak lines 
on the linear part) as well as increasing strengths (i.e., extremely strong damping-dominated 
lines in the damping part), reflecting the importance of the core-deepening effect caused by the 
dilution of 5'l as already mentioned in subsection 2.2 (cf. figure 2). 

This naturally explains the behaviors of the extent of the non-LTE corrections for each line 
at different metallicities; for example, the reason why the weak 6154/6161 lines in disk stars 
show almost negligible non-LTE effects (cf. figure 5c), why the |A| values for the 8183/8195 
attain maximum at [Fe/H] ^ — 1 (cf. figures 5d, 6b, and 7b), and why the |A| values for the 
5890/5896 lines systematically decreases from [Fc/H] ^ —2 to ^ —4 (cf. figure 7b). 

We should note, however, that the situation more or less differs from line to line, refiecting 
their individual properties. For example, lines at longer wavelengths tend to suffer larger non- 
LTE effects because of the decreasing sensitivity of the Planck function to the temperature 
(i.e., the minimum allowed LTE residual intensity tends to be raised and thus the difference 
between LTE and non-LTE becomes more significant). Hence, we can understand the fact that 
the near-IR 8183/8195 lines show relatively larger non-LTE effects compared to the other lines 
at similar line-strengths. 
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Figure 6 (a) [Na/Fe] vs. [Fe/H] relation resulting from the non-LTE re- 
analysis of the rather old data published before 1990. See panel (b) for the 
data sources of equivalent widths and the corresponding symbols. Lines show 
the representative theoretical predictions: Dash-dotted line (S98) taken 
from figure 8 of Samland (1998), dashed line (TWW) — taken from figure 
17 of Timmes et al. (1995) corresponding to the two cases for the standard 
Woosley and Weaver's (1995) Fe yield from massive stars and the reduced 
one by a factor of 2 (which may be more reasonable according to their sug- 
gestion), dashed line — taken from figure 7 of Goswami and Prantzos (2000), 
for the two cases of Na yield; time-independent constant yield [GP (const); 
only for a comparison purpose] and the realistic metallicity-dependent Na 
yield [GP(var)]. (b) The corresponding NLTE corrections used for deriving 
the [Na/Fe] values shown in panel (a), plotted as functions of [Fe/H]. 
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Figure 7 (a) [Na/Fe] vs. [Fe/H] relation resulting from the non-LTE re- 
analysis based on the recently published new data after 1990. The data point 
for CS 22949-037 ([Na/Fe] = -hl.94 at [Fe/H] = -3.99) derived from the re- 
analysis of McWilliam et al.'s (1995a) EW data is outside the plot range 
of this figure (cf. figure 8a for a more global view including this anoma- 
lous point), (b) The corresponding NLTE corrections used for deriving the 
[Na/Fe] values shown in panel (a), plotted as functions of [Fe/H]. Otherwise, 
the same as in figure 6; see the caption therein for more details. The re- 
sults of the reanalysis of Pilachowski et al.'s (1996) 5890/5896 data for the 
low-gravity giants/supergiants (filled triangles) may suffer rather large un- 
certainties and should not be seriously taken (cf. the footnote in subsection 
5.3). 
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Figure 8 (a) [Na/Fe] vs. [Fe/H] relation resulting from the non-LTE re- 
analysis of the new literature data after 1990. The plotted data are essentially 
the same as in figure 7a, but dwarfs and giants are discriminated by different 
symbols: open circles — high-gravity stars (log g > 3), filled circles low- 
gravity stars (logg < 3). The theoretical prediction of Timmes et al. (1995) 
(the case of reduced Fe yield denoted by "TWW (1/2 x)"; cf. the caption of 
figure 6a), which appears to match the observed trend most satisfactorily, is 
also drawn for a comparison, (b) Tgs vs. log 5 relation for the sample stars 
relevant to the data plots in (a). The meanings of the symbols as the same 
as in (a). 
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Based on figures 5-7, we can assess the importance of the non-LTE effect on different hnes 
in sodium abundance determinations as follows: The 6154/6161 lines are the best abundance 
indicators, for which the non-LTE effect is practically insignificant (i.e., ;^ 0.1 dex) and may 
be neglected, though they are usable/visible only at —2 ^ [Fe/H] ^ 0.5). The 5683/5688 lines 
may also be tolerably analyzed with LTE at the metallicity of —3 ^ [Fe/H] ^ —1. Regarding 
the 5890/5896 and 8183/8195 lines, which are important for diagnosing the Na abundance 
behaviors of very metal-deficient stars because of their strengths (especially, the 5890/5896 D 
lines are the only tool for investigating the extremely metal-poor regime of —4 ^ [Fe/H] ^ —3), 
we consider that properly taking account of the non-LTE effect is necessary. 

5.2 Disk Stars 

The overall behavior of the Na abundances for disk stars (— 1 ^ [Fe/H] ^ +0.4) has long been 
known to nearly scale with Fe; i.e., [Na/Fe] ~ (e.g., Wallerstein 1962; Tomkin et al. 1985). 
There remains almost no doubt for this broad description, which may be confirmed also in our 
figures 5-7.'^^ 

However, our attention is paid rather to the more delicate systematic trend in the [Na/Fe] 
vs. [Fe/H] relation; namely, the "upturn" revealed by Edvardsson et al.'s (1993) analyses based 
on the 6154/6161 lines, characterized by slight upward increases on both sides of the broad 
minimum around [Fe/H] ~ —0.2. Carretta et al.'s (2000) reanalysis on Edvardsson et al.'s 
(1993) data also reproduced this tendency, as shown in figure 9 in their paper. Similarly, our 
reanalyzed results on Edvardsson et al.'s (1993) data are essentially the same as their original 
results (cf. small open circles in figure 5c). 

Supportive evidences for this trend have been reported from successive analyses on in- 
dependent observational materials. Feltzing and Gustafsson (1998) confirmed in their Na i 
6154/6161 analyses of metal-rich stars that [Na/Fe] progressively increases from ~ ([Fe/H] 
0) to +0.2 ([Fe/H] ~ +0.4) (see the dots in figure 5c). They also investigated its behavior 
in the kinematically different stellar groups to see if there is any position-dependent effect in 
the Galactic chemical evolution, though distinct differences were not observed. Also, Chen et 
al.'s (2000) LTE analyses of 6154/6161 lines based on the BAO data of disk stars (common to 



ttNotc that figures 6-8 include low-gravity giants of population I [e.g., reanalysis results of Korochin and Mishen- 
ina's (1999) data shown by open stars in figure 7] showing markedly positive [Na/Fc] values in contrast to 
others. These stars arc expected to have suffered evolution-induced Na cnrichnient in the envelope caused by 
the drcdge-up of NeNa cycle products in the H-burning shell, such as the case of supergiants (Takeda, Takada- 
Hidai 1994). Hence, their abundance characteristics should be discriminated from those related to the Galactic 
chemical evolution. 
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this study) suggested a weak minimum at [Fe/H] —0.2 superposed on the nearly flat [Na/Fe] 
(~ 0) (cf. figure 8 therein) consistent with that found by Edvardsson et al. (1993). 

As mentioned in the previous subsection, the non-LTE corrections on the Na i 6154/6161 
lines usually used in sodium abundance analyses of disk stars are so small that they are prac- 
tically negligible. Therefore, it is expected that our non-LTE analyses on these orange lines 
of multiplet 5 reproduce essentially the same results as obtained by Chen et al. (2000), which 
is actually confirmed in figure 5c. We point out, however, that our non-LTE [Na/Fe] values 
(on the BAG data) derived from other stronger lines (5683/5688, 5890/5896, and 8195) also 
show signs of "upturn", (tough not so clear as the case of 6154/6161 lines) in spite of their 
inadequacy for abundance determinations, as can be seen from figures 7a, b, and d. Hence, we 
may state that this trend is firmly ascertained. 

Considering the decline of the [Na/Fe] values for a decrease of [Fe/H] both on the metal- 
poor side ([Fe/H] <. —1; cf. the next subsection) and the metal-rich side (0 <. [Fe/H]; cf. 
Feltzing, Gustafsson 1998), the existence of such an "upturn" may indicate an extra production 
of sodium at the cosmic metallicity of — 1 [Fe/H] ^ —0.5. However, any speculation had 
better be put off until the behaviors of other elements closely related to the synthesis of Na, 
especially Al and Mg, have been established. 

According to our reanalysis results for Prochaska et al.'s (2000) data (cf. diamonds in figure 
7), the [Na/Fe] values of thick disk stars do not appear to exhibit significantly distinct differences 
from those of normal thin-disk stars, though our mean non-LTE ([Na/Fe]) averaged over their 
10 stars are nearly zero [-0.03(a- = 0.06) and -|-0.05(£7 = 0.04) for 5683/5688 and 6154/6166 
lines, respectively] and slightly smaller than the value (0.087) they obtained, which may be due 
to the non-LTE effect (mean non-LTE corrections are —0.09 and —0.06, respectively). 

5.3 Halo Stars 

Now, wc discuss the behavior of [Na/Fc] for mctal-poor halo stars. For simplicity, our discussion 
is confined to figure 7, the results of our reanalyzes based on the data taken from the works 
after 1990, which may presumably be more accurate and reliable than those in figure 6 for the 
reasons mentioned in subsection 4.3. 

One important feature that can be read from figure 7 is the systematic decrease of [Na/Fe]; 
i.e., from [Na/Fe] '--^ at [Fe/H] ~ —1 down to [Na/Fe] ~ —0.4 (though a rather large diversity 
of ~ ±0.2) at [Fe/H] ^ —2. This means that we have reasonably reconfirmed the recent results 
of Baumiiller et al. (1998) and Stephens (1999), as opposed to the previous belief imtil mid 
1990's (cf. section 1). It should be stressed that each of the different multiplet lines [5683/5688 
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(Stephens et al. 1999; half-filled circles), 5890/5896 (McWiUiam et al. 1995a; filled circles), 
6154/6161 (Kraft et al. 1992; double circles) and 8183/8195 (our OAO data; St. Andrew's 
crosses)] yield consistent results with each other, in spite of the considerably different extents 
of the non-LTE corrections. Thus, we may state that this subsolar [Na/Fe] ratios in metal-poor 
stars of —3 ^ [Fe/H] ^ —1 are firmly established,** while previous LTE analyses using the 
non-LTE sensitive lines (e.g., McWilliam et al. 1995b, who invoked 5890/5896 lines) must have 
overestimated the sodium abundances. 

Turning our attention to the further extreme metal-poor regime of —4 ^ [Fe/H] ^ —3, 
we see in figure 7 an interesting trend of rising [Na/Fe] for a lowering of the metallicity, again 
recovering [Na/Fe] ^ at [Fe/H] ~ —4, though a rather insufficient number of stars in this 
region prevents us from making any definitive argument. It is highly desired to increase the 
sample of those stars toward establishing the behavior of [Na/Fe] in this ultra-low metallicity 
region, which is very important for investigating the history of the early-time Galaxy as well as 
for constructing the realistic model of Galactic chemical evolution. 

Finally, we compare the observed [Na/Fe] vs. [Fe/H] relations obtained for the metal- 
deficient halo stars with representative theoretical predictions calculated by using various chem- 
ical evolution models of our Galaxy. As mentioned in section 1, since the production of Na is 
dependent on the excess neutrons, its yield should naturally be metallicity-dependent and in- 
crease with time. As a matter of fact, the assumption of constant yield results in a totally 
unrealistic [Na/Fe] vs. [Fe/H] relation with a markedly negative rf[Na/Fe]/d[Fe/H] gradient 
(cf. the dashed line labeled "GP(const)" in figures 6 and 7, which was taken from figure 7 of 
Goswami, Prantzos 2000). It may be possible, therefore, to judge which modeling of Na produc- 
tion (among several proposed chemical evolution calculations) represents the actual situation 



•'■'■Some remark may be due regarding our non-LTE rcanalysis results of Pilachowski et al.'s (1996) 5890/5896 data 
(cf. filled triangles in figure 7) for population II lo'w-gravity giants mostly having parameter values of T^b ~ 
4000-5000 K, logs ~ 0-2, and —2.5 < [Fe/H] < —1.5, which do not reveal any such clear subsolar tendency 

[([Na/Fe]) = 0.00(±0.24) on the average] unlike the other cases mentioned above, despite their original LTE 
results showed a moderate subsolar trend (([Na/Fe]) = —0.17; cf. Pilachowski et al. 1996). We point out 
that some of their program stars are low temperature/gravity K giants with parameters of Teg ~ 4000 K or 
logg ~ 0.0, which arc outside the parameter grid of our non-LTE correction tables and the applied non-LTE 
correction had to be evaluated by an extrapolation. In these low temperature/gravity cases, it exceptionally 
happens that the uou-LTE effect tends to act as a line- weakening mechanism (i.e., marginally positive non-LTE 
correction) due to the enhanced Si^ {Si^/B > 1) in the line- forming region according to our calculations, as 
mentioned in footnote 1 in subsection 2.2. Hence, though appreciable positive non-LTE corrections are actually 
observed for several stars in figure 7b (filled triangles), these are all such low-gravity K giants for which this 
effect may have been rather erroneously exaggerated by inaccurate extrapolations. Also, because of the low- 
temperature nature of the Pilachowski et al.'s (1996) program stars, the strengths of the 5890/5896 D lines are 
considerably large (200-400 mA) in spite the low metallicity and comparatively large ambiguities arc expected 
(e.g., equivalent-width measurements, how to choose the damping parameter or the microturbulent velocity). 
Consequently, our reanalysis results on their data may suffer large uncertainties and should not be seriously 
taken. 
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best by comparing the predicted [Na/Fe] vs. [Fe/H] relation with that observed for halo stars 
mentioned above. 

For this purpose, theoretical predictions taken from three recent representative works are 
depicted in figures 6 and 7; i.e., those of Timmes et al. (1995; solid lines), Samland (1998; 
dashed-dotted line), and Goswami and Prantzos (2000; dashed lines). We can draw from figure 
7 the following conclusions. 

— Among these, Timmes et al.'s (1995) results, especially for the case of reduced Fe yield 
denoted as "TWW(l/2 x)" (cf. figure 8a), appear to most satisfactorily reproduce the observed 
tendency of [Na/Fe] (i.e., broad/shallow dip- like feature with a minimum around [Fe/H] ~ —2) 
mentioned above. Their calculations are based on Woosley and Weaver's (1995) yields, simple 
dynamical model for the Galaxy with infall, a standard Salpeter (1955) initial-mass function 
(IMF) and a quadratic Schmidt star-formation rate. (Each line corresponds to the different 
choice of the Fe yield from massive stars; i.e., that from the original Woosley and Weaver paper 
and that reduced by a factor of 2; the latter is rather recommended by Timmes et al. 1995). 

— Samland's (1998) model is based also on Woosley and Weaver's (1995) yields. However, 
since he adjusted the metallicity-dependent Na yield in his own way [cf. equation (9) therein] 
so as to fit the observed trend of [Na/Fe] known at the time of mid-1990's (i.e., erroneously 
overestimated to be [Na/Fe] ~ even for halo stars), his prediction does not fit the observed 
[Na/Fe] obtained in this study. 

— Goswami and Prantzos's (2000) calculation (labeled "GP(var)"in figure 7) explains the 
decline of [Na/Fe] from [Fe/H]~ —1 to [Fe/H]~ —2 reasonably well. However, their prediction 
suggests an ever-decreasing (or a plateau-like) [Na/Fe] toward the extremely low metallicity 
regime ( [Na/Fe] ~ —4), which apparently contradicts the tendency we found. While their model 
is based on Woosley and Weaver's (1995) metallicity-dependent yields (as was done by Timmes 
et al. 1995) with a Fe yield reduced by a factor of 2, they used a different IMF (Kroupa et 
al. 1993) and a different halo model. From this point of view, the dynamical model of the 
Galaxy and IMF adopted by Timmes et al. (1995) might be preferable to those of Goswami 
and Prantzos (2000), as far as the early Galaxy ([Fe/H] ^ —2) is concerned. 

6 CONCLUSION 

Wc carried out extensive non-LTE calculations on the atomic model of neutral sodium and 
model atmospheres with a wide range of parameters, for the purpose of non-LTE abundance 
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determinations based on eight representative Na i lines at 5683, 5688, 5890, 5896, 6154, 6161, 
8183, and 8195 A. 

The non-LTE effect almost always acts as a line-strengthening mechanism (thus the non- 
LTE abundance correction being negative). Though its extent differs from line to line (from 
almost negligible level of a few hundredths dex to a considerable amount of ~ 0.5 dex), strong 
saturated lines (though not too strong to be damping-dominated) tend to suffer a large non- 
LTE effect. Generally speaking, the lines for which non-LTE corrections had better be taken 
into account in Na abundance determinations are 5890/5896 and 8183/8195 lines, while the 
5683/5688 and (especially) 6154/6161 lines are comparatively less sensitive to a non-LTE effect. 

Our non-LTE reanalyzes of our BAO and OAO equivalent-width data along with those 
taken from extensive literature revealed the following conclusions: 

— (1) Regarding the disk stars with —1 ^ [Fe/H] ^ 0.4, we confirmed the existence of a delicate 
"upturn" feature (i.e., broad/shallow dip around the minimum at [Fe/H] ~ —0.2) superposed 
on the general tendency of [Na/Fe] ~ 0, as has been reported by recent analyses on weak Na i 
6154/6161 lines being inert to any non-LTE effect. Moreover, according to our analyses on our 
BAO data, even the abundances derived from strong (saturated or damping-dominated) lines, 
which are generally unsuitable for abundance analyses, suggested this tendency. 

— (2) We found based on our analyses of recent data after 1990's that the [Na/Fe] ratios 
for metal-poor halo stars show a "subsolar" behavior; i.e., [Na/Fe] decreases from ~ (at 
[Fe/H]'~ —1) to ~ —0.4 (at [Fe/H]~ —2), while it appears to rise again with a further decrease 
in metallicity toward recovering [Na/Fe]'~ again at [Fe/H]~ —4. It is evident that the previ- 
ous suggestion of almost solar Na-to-Fe ratio over all metallicity, which was actually believed 
until mid 1990's, is attributed to the overestimation of sodium abundances due to the neglect of 
non-LTE effects. Among the representative theoretical models of Galactic chemical evolution, 
that of Timmes et al. (1995) appears to reasonably reproduce this behavior. 

— (3) Given this observational evidence, the run of [Na/Fe] is in rough accord with that of 
[Al/Fe] (cf. Baumiiller, Gehren 1997), which would be a sound consequence from a theoretical 
point of view, since both are expected to show similar behaviors (see, e.g., subsection 5.8 of 
Samland 1998). In order to further clariiy the connection between the present results and the 
abundances of such important Na-related elements in metal-poor stars, we are now planning to 
carry out extensive non-LTE analyses on Al and Mg, in the way similar to that adopted in this 
paper by using our own as well as the literature data, the results of which will be reported in 
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our forthcoming papers. 
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